Anderson impurity in the one-dimensional Hubbard model on finite size systems 
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An Anderson impurity in a Hubbard model on chains with finite length is studied using the 
density-matrix renormalization group (DMRG) technique. In the first place, we analyzed how the 
reduction of electron density from half-filling to quarter-filling affects the Kondo resonance in the 
limit of Hubbard repulsion [7 = 0. In general, a weak dependence with the electron density was 
found for the local density of states (LDOS) at the impurity except when the impurity, at half- 
filling, is close to a mixed valence regime. Next, in the central part of this paper, we studied the 
effects of finite Hubbard interaction on the chain at quarter-filling. Our main result is that this 
interaction drives the impurity into a more defined Kondo regime although accompanied in most 
cases by a reduction of the spectral weight of the impurity LDOS. Again, for the impurity in the 
mixed valence regime, we observed an interesting nonmonotonic behavior. We also concluded that 
the conductance, computed for a small finite bias applied to the leads, follows the behavior of the 
impurity LDOS, as in the case of non-interacting chains. Finally, we analyzed how the Hubbard 
interaction and the finite chain length affect the spin compensation cloud both at zero and at finite 
temperature, in this case using quantum Monte Carlo techniques. 

PACS numbers: 73.23.-b, 73. 21. La, 72.15.Qm, 71.27.-|-a 



I. INTRODUCTION 

There is an increasing interest in strongly correlated 
systems containing onedimensional (ID) structures, both 
for their novel properties as well as for their possible ap- 
plications. For example, chains and two-leg ladders are 
present in transition metal oxides and in organic con- 
ductors. Other systems of increasing interest involv- 
ing strongly correlated features are nanoscopic devices 
which have been recently fabricated and experimentally 
studiedi In many of these devices, a droplet of electrons 
confined in the three spatial directions, the quantum dot 
(QD), is connected to metallic leads. In this case the 
leads are described by tight-binding models, or in gen- 
eral by Fermi liquids. The transport through such de- 
vices involve the formation of a Kondo resonance in the 
QD. Kondo physics is then relevant to understand this 
kind of devices r~ 

It is also possible to connect the QD to conducting 
leads with electron correlations. This is the case of spin 
valves, one of the simplest devices used in spintronicS)2i4 
where a QD is attached to ferromagnetic leads. These 
leads can be built from materials such as manganese 
oxides^ or other compounds containing transition metal 
oxides such as vanadatesi^ There have also been recent 
developments in producing devices where the leads are 
formed by carbon nanotubes (although ultimately these 
leads are connected to metallic terminals) which can be 
considered quasi-onedimensional and where electron cor- 
relations are presumably important iSiS In fact, a Lut- 
tinger liquid behavior, typical of correlated ID systems, 
has been found in carbon nanotubesfifi 

The problem of a magnetic impurity or QD in a cor- 
related onedimensional system is then an interesting 
topic both for fundamental and for technological reasons. 



From the theoretical point of view, the most important 
results for this kind of systems were obtained by con- 
sidering a single impurity attached to leads described by 
Luttinger liquids. Much less is known when the leads are 
described by a Hubbard model. In this case, this model 
has been treated by linearizing the dispersion relation 
and, using bosonization techniques, the long wavelength 
behavior is captured. Within this approach, it was shown 
that in a system of a magnetic impurity connected to 
Luttinger chains by a Kondo exchange interaction J, the 
Kondo temperature acquires the power-law expressionii 
Tk ~ ( J)2/(i--'^p) instead of the conventional exponen- 
tial law for noninteracting chains. The Luttinger expo- 
nent Kp also determines the conductance through an in- 
teracting wire, which is given by G = 2K pe^ / 

In this article we consider the even less studied prob- 
lem of a single-impurity Anderson-Hubbard Hamiltonian 
on finite-size systems. Our aim is to provide a detailed 
microscopic picture of the competition between the elec- 
tron correlations on the impurity and the ones on the 
leads, and to understand the role of finite length of the 
chains. Although the lengths considered in this study are 
still smaller than the physical length of current devices, 
the study on finite-size systems are interesting per se. In 
fact, it has been suggested that the Kondo temperature 
can be very different in finite chains with respect to the 
one in infinite leadSfiS, It is also expected that the size of 
nanoelectronic devices will be further reduced in the near 
future. The fundamental problem of the competition be- 
tween the Kondo effect and the Luttinger liquid starting 
from a microscopic Hamiltonian was recently studied on 
infinite systems using an approximated analytical renor- 
malization group technique^ The problem of an Ander- 
son impurity in a t — J model has been also addressed 
within the Bethe Ansatz formalism»i^ 

Although manganese oxides should be described by 



2 



a generalized ferromagnetic Kondo-lattice modet--, and 
vanadates or carbon nanotubes by some variants of the 
Hubbard or i — J models on laddersi^iiS, in order to 
get an insight on the physical behavior, it is necessary 
to consider first the simpler case of leads described by 
a ID Hubbard model. Of course, some of the results 
obtained are relevant to magnetic impurities in strongly 
correlated quasi-lD compounds such as those mentioned 
at the beginning. Results for LDOS for example, could 
be compared with those obtained by scanning tunneling 
microscopy (STM) experiments.-^^ 

Numerical techniques are ideally suited to cope with fi- 
nite systems. In particular we mainly employ the DMRG 
technique^ which provides essentially exact results for 
various real-space properties on finite chains. These real- 
space properties of real-space models can shed light on 
the functioning of QD devicest^ In this sense, in this 
work we will determine the way in which the Kondo ef- 
fect is affected by varying the electron density or the 
Hubbard repulsion on the leads by computing the elec- 
tron occupation and the square of the ^-component of 
the spin at the impurity site. More importantly, from 
the spin-spin correlation functions we will study the "spin 
compensation cloud" , a possible measure of the more elu- 
sive Kondo screening cloud. Using this technique, it 
has been shown that the screening cloud is reduced by 
electron correlations on the leads in a system with two 
Kondo impurities. ^'^ We will extend this approach to the 
Anderson impurity where we will make a careful study 
of finite-size effects. DMRG results for the spin compen- 
sation cloud have been reproduced to a good approxima- 
tion by using a recently developed quantum Monte Carlo 
technique^. This technique also allows us to study this 
quantity at finite temperature. 

In order to study transport properties, we will study 
the LDOS at the impurity, which is related to the conduc- 
tance in an essential way. By implementing some recent 
developments in DMRGSL2&, we will compute the con- 
ductance as the response of the system to a finite but 
small voltage bias applied to the leads. We will examine 
this property as the Hubbard repulsion U is varied. We 
will also show that the effect of this variation on both the 
LDOS at the impurity and the conductance, depend in 
turn on the values of the electron interactions on the im- 
purity. We will explore various types of impurities, from 
the mixed valence to Kondo regimes. 

The paper is organized as follows. In Section ^ we 
describe the model studied and we provide details of the 
numerical techniques employed. In Section IlIII we study 
the effect of reducing the electron density from half-filling 
to quarter-filling for the case of noninteracting (tight- 
binding) leads. In Scction llVl we analyze how the Kondo 
effect varies with the Hubbard repulsion U on the leads at 
quarter-filling. Finally, the relation between the results 
obtained and those in the previous literature, and the rel- 
evance of the present study to real systems or nanoscopic 
devices, are discussed in Section Ivl 



II. MODEL AND METHODS 

In this paper we use the terms "impurity" and QD 
indistinctly, and we refer to the remainder sites of the 
chain as the "leads". We consider a one-dimensional 
single-impurity Anderson-Hubbard model defined by the 
Hamiltonian: 

n = - t ^ {cl^c^+i„ + H.c.) + U ^ "j,T»^<,i 

i=<-2,a i=<-l 

~ * X! i4aCt+la + H.C.) + "-»,T'^«a 
i=>l.,<y j=>l 

- t' ^{c^-laCGo + cl^Cia + H.C.) 
a 

+ e' no -I- C/'7iojno,x (1) 

where conventional notation was used. The first two 
terms correspond to the left lead and the following two 
terms to the right lead. The next term is the impurity- 
lead interaction, and the last two terms are the on-site 
energy and Hubbard repulsion at the impurity site, lo- 
cated at site 0, respectively. We adopt t as the scale of 
energy. In the case oiU — 0, the Hamiltonian J^l reduces 
to the single-impurity Anderson model, and for d = 0, 
[/' = U,t' ^ t, to the Hubbard model. 

Hamiltonian ^ will be studied mostly using the nu- 
merical technique DMRGS2*2^ on finite-size clusters of 
length L with open boundary conditions (OBC). This 
algorithm provides numerically exact results for static 
properties at zero temperature with a precision which de- 
pends on the number M of states retained. Most of the 
results here reported were obtained for M = 600, except 
otherwise stated, asserting us that the integrated weight 
of discarded states are of order 10~^ in the worst case. 
On the other hand, results for dynamical properties, such 
as LDOS, should be taken qualitatively, as discussed be- 
low. The same apply to the results obtained for the time 
evolution of the current on the links connecting the im- 
purity with the leads. In order to assess even-odd effects 
most of the results reported below were obtained by using 
even and odd chain lengths. Throughout all this paper, 
DMRG calculations will be done in the subspace of total 
5;, = {S^ = 1/2) for L even (odd). The QD is located 
in one of two central sites of the chain when L is even (in 
the central site when L is odd). 

We have computed static properties such as the elec- 
tron occupancy on each site, {fii^a) (f =Tii) and spin- 
spin correlations from the QD, S{j) = {S^ SI)-{S^){S'-). 

A very important quantity related to the Kondo effect 
is the Kondo screening length, which is somewhat elusive 
to computeJ^ A possible measure of this Kondo length 
is the length of the "spin compensation cloud"— defined 
as the length ^ such that 

S{j) = xS{0), (2) 
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where x is an arbitrary parameter. In the following we 
adopt a fixed value x — 0.9 in order to compare ^ for the 
different cases. 

One of the main quantities we have studied is the lo- 
cal density of states, p{uj), at the QD. In the first place, 
from this quantity it would be possible to evaluate the 
conductance in the linear response regime (see below). 
In the second place, recent advances in STM have made 
it possible to directly measure this quantity. Since in 
our DMRG calculation the measurements are performed 
when the two added sites are at the center of the chain 
(symmetrical configuration) , the QD then is one of these 
two sites which are exactly treated. Then, we adopt the 
approximation of applying the creation and annihilation 
operators at the QD on the ground state vector and then 
determine p{w) following the well-known continued frac- 
tion formalism. A more accurate approach would be, 
after the application of each of those creation and an- 
nihilation operators, to run additional sweeps for an en- 
larged density matrixi2& In any case, the truncation of the 
Hilbert space is the essential source of error in DMRG, 
and to estimate the precision of our approach we have 
compared results for various numbers of retained states 
M, from M = 300 to 600. We would also like to stress 
the fact that the conductance in linear response is re- 
lated to the LDOS near w = where the approximation 
is more precise. Our calculation starts to be less precise 
as we move away from the chemical potential, since then 
higher excited states are involved. 

For this kind of systems, specially in connection to 
nanoscopic devices, the most interesting properties are 
those related to transport, in particular, the conductance 
G. When the chains or leads are described by a tight- 
binding model, or in general by a Fermi liquid, various 
expressions relating it to the LDOS at the impurity are 
available in linear response, as mentioned before. In ad- 
dition, the Friedel's sum rule (see Eq. iQ below) allows a 
simple and precise calculation of G, particularly for nu- 
merical techniques. In the case of Hubbard chains, or 
in general Luttinger liquid leads, one must resort to the 
expression, in linear response^i^ 

G = lim / [ dxdt e''^'{J{x,t)J{0,0)) (3) 

uj^o HLuj Jq J 

where the current is defined as 

J{x,t) = t^t, E^^WK'^l+i-^- - H.c.)\^{t)), (4) 

and |^(i)) is the ground state of the system at time t. 

Although such a calculation can be carried out within 
DMRG, some recent developments collectively known as 
time-dependent DMRG2L2& allow to compute directly 
{J{x,t)) as {'^{t)\J{x)\'^{t)). The time-evolution is trig- 
gered by the application of a bias voltage AV between 
the left and right leads a,t t = OiSSiffi That is, at time 
t = 0, a potential —AV/2ni {AV/2ni) on the sites on 
the left (right) leads is switched on. Then, charge moves 



from one lead to the other in a non-equilibrium process. 
This non-equilibrium process was studied analytically in 
Ref. 1^ In finite size systems, {J{x, t)) follows an oscilla- 
tory evolution as charge is bounced between the left and 
right leads. As the length of the leads are extended to 
infinity, the period of this oscillation becomes infinitely 
long (the leads become infinite charge reservoirs) and the 
usual co ncept o f conductance is recovered. It was shown 
in Refs. I29l30l that even for moderate size clusters, on 
each period, {J{x,t)) reaches a plateau that corresponds 
to the value of the bulk limit. Hence a measure of the 
conductance can be obtained as the value of this plateau 
or the amplitude of this oscillation, that is: 

G(AF) = max^ (5) 

where 

J{t) = {JL{t) + JR{t))/2 (6) 

and Jhit) (JRit)) is the current on the link connecting 
the impurity site with the left (right) lead. Notice that 
Eq. ((Sj is valid for an arbitrary finite bias AV. The evo- 
lution of the ground state |^'(i)) from t to t + At is com- 
puted with a static Runge-Kutta treatment of the time- 
dependent Schrodinger equationi^ Again in this case, a 
much precise calculation could be performed by includ- 
ing |^(i)) at intermediate times r as targeted states in 
the density matrix and by running additional sweeps 
This approach is of course much more expensive compu- 
tationally. 

In the calculations reported below we adopted a bias 
AV = 0.01 and a time step At = 0.1. 

Finally, we report finite temperature results ob- 
tained with quantum Monte Carlo simulations using the 
stochastic series expansion (SSE) algorithm. In this 
case, we consider chains with periodic boundary condi- 
tions (PBC). The use of PBC helps to assess the effects 
of the Friedel oscillations which appear with OBC, used 
with DMRG, but are absent with PBC. SSE with the loop 
upgrade is a very efficient QMC algorithm which allows 
to reach much lower temperatures and larger sizes than 
other QMC algorithms. It also avoids the problem of ex- 
trapolation of the Trotter dimension common in "world- 
line" algorithms. It is also worth to note that SSE works 
in the grand canonical ensemble, i.e., the chemical po- 
tential has to be tuned to give a desired density on aver- 
age. Notice that the "minus sign problem" which affects 
quantum Monte Carlo simulations for fermionic models 
is absent in ID with hopping between nearest neighbor 
sites only since fermion permutations appear just as a 
boundary effect. 

III. NONINTERACTING LEADS, ELECTRON 
DENSITIES AWAY FROM HALF-FILLING 

The main purpose of this work is to study the effects of 
correlations on the leads. This study will be performed 
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FIG. 1: (Color online) (a) Electron occupancy and (b) (5,2) 
at the impurity site as a function of the filling for U' = 8, 
t' = 0.5 (circles), U' = 8, t' = 1.0 (squares), U' = 4, t' = 0.5 
(diamonds) and U' = 2, t' = 0.5 (triangles). L = 96 



in the next Section. Most theoretical studies on Kondo 
effect have considered the half- filled case [n — 1), where 
the system is particle-hole symmetric. Now, since an in- 
finitesimal value of U drives the system into an insulating 
state, in order to keep the leads metallic in the presence 
of correlations it is necessary to study the system away 
from half- filling. In Section llVI we will work at quarter- 
filling {n = 0.5). In this Section, as a preliminary step, 
we study the evolution of the Kondo effect as the filling is 
reduced from n = 1 to n = 0.5 in the absence of interac- 
tions on the lead ([/ = 0). In particular we arc interested 
in following the evolution of what at half-filling is the 
Kondo resonance, as the electron density decreases. 

Model Eq.Q, which as stated above reduces to the 
single-impurity Anderson model, was studied by DMRG 
for several L up to 96 sites, and for fillings n — \, 0.875, 
0.75, 0.625 and 0.5. Since, as we said in the Introduc- 
tion, there is a competition between the correlations on 
the leads and the correlations at the impurity or QD, 
we consider four sets of interactions on the impurity: 
{{/' = 8,t' = 0.5}, {[/' = 8,<' = 1.0}, {[/' = 4,i' = 0.5}, 
and {U' ~ 2,t' = 0.5}. In all cases we consider the 
symmetric case, e' — —U'/2. For these sets of parame- 
ters, the effective Kondo coupling Jeff = it'^/U' takes 
the values 0.125, 0.5, 0.25, and 0.5 respectively, although 
strictly this relationship is only valid for U' >> t' , and 
hence not much applicable to the last case. In spite of 
sharing the same Jeff, the second and fourth sets behave 
quite differently as we will show below. 

Let us start by examining the dependence of the elec- 
tron occupancy (no) and {S'^ q) at the impurity site with 
the electron density. This variation depends in turn 
strongly on the interactions at the impurity as it can 
be seen in Fig.^ If (no) = 1, {S'^g) takes its maximum 
value (1/4) when the impurity is a perfect spin-1/2, corre- 
sponding to the full Kondo regime. In the mixed valence 
regime, it takes the value 1/12. For the parameter set 
{U' = 8,t' = 0.5}, the system is in a well defined Kondo 
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FIG. 2: (Color online) (a) Conductance obtained from the 
Friedel's sum rule, and (b) length of the compensation cloud 
^, as a function of the filling for U' = 8, t' — 0.5 (circles), 
U' = 8, = 1.0 (squares), U' ^ 4, t' = 0.5 (diamonds) and 
U' = 2, t' = 0.5 (triangles). L = 64 (full symbols, dot lines), 
L = 96 (open symbols, full lines). 

regime, and the dependence with n is relatively weak. 
On the opposite case, for {[/' = 2,t' = 0.5} the system 
is not in a well-defined Kondo regime at half-filling, and 
the dependence with the electron filling is stronger. The 
reduction of electron density drives the system into the 
mixed valence regime. Notice that in this second case, 
the occupancy of the impurity is smaller than n. Besides, 
it is well known that there are Friedel oscillations in (rii) 
beginning at the ends of the chains with OBC.'^^ It is then 
valid to ask if these oscillations could have some effect on 
the electron occupancy at the impurity. In the present 
case, any important influence of the Friedel oscillations 
on the impurity properties can be ruled out because the 
results for L = 96 are virtually identical within error 
bars to those obtained with SSE and using PEC where 
these oscillations are absent, at the lowest temperatures 
we considered. Of course, the presence of the impurity 
induces other density oscillations which overimpose to 
the Friedel oscillations induced by the open ends. In fact 
we have observed in some cases that the presence of the 
impurity changes the period of the open end oscillations 
from 2k F to 4k p. This issue is worth to be examined 
more thoroughly but it is somewhat outside the scope of 
the present work. 

From the values of the electron occupancy it is possi- 
ble to obtain the conductance by using the Friedel's sum 

Ga — sin[-Kna)'^ (7) 

(cr =}?!) which is valid when the bandwidth is much 
larger than the Kondo coupling. Results for G = Gj + G| 
are shown in Fig. |2Ia), for the four types of impurities 
considered. As expected from the results in Fig.^a), the 
conductance decreases by reducing electron density from 
half-filling. 

The length of the spin compensation cloud ^ is de- 



picted in Fig.[5Jb) for the same four types of impurities, 
f depends on two factors: (i) the effective exchange cou- 
phng between the magnetic impurity and the spins of 
the conduction electrons, and (ii), to a lesser extent, how 
much defined is the magnetic character of the impurity 
or the Kondo regime. The finite length of the cluster im- 
poses an upper bound to the compensation length. It has 
been pointed out that the inability of the finite system 
to accommodate the Kondo cloud would modify in turn 
the Kondo effect, in particular the Kondo temperatureii^ 
For the parameter set U' = 8, t' ^ 0.5 (J^// = 0.125) 
the compensation cloud is large, quite likely exceeding 
the system size since ^ w L/2, and scaling linearly with 
the system size. For U' = 4, = 0.5 {Jeff = 0.25), ^ 
is somewhat smaller but still scales linearly with L. On 
the other hand, for the U' = 8, t' = 1 {Jeff = 0.5) im- 
purity, where {S'^ q) « 0.21, the length of the compensa- 
tion cloud is much smaller and it has a weak dependence 
with L. An intermediate situation occurs for the fourth 
type of impurity, U' = 2, t' = 0.5 which has the same 
Jeff = 0.5 but with a poorly defined magnetic character 
{{Slo) < 0.2), leading to a more extended compensation 
cloud. Notice an overall weak dependence of ^ with the 
electron density. 

It is important to understand the behavior of the con- 
ductance shown in Fig.E^a) by studying the LDOS at the 
impurity site, p{uj). In Fig. 13 we show p{lu) close to the 
Fermi level, for the four types of impurities and for several 
electron densities from half-filling to quarter-filling. p{u!) 
has been shifted to the chemical potential p for the sake 
of comparison, since starting from half-filling, fj. shifts to 
Lu < 0. In this figure and all the following similar ones, we 
adopted a Lorentzian broadening of the peaks of S — 0.1. 
A measure of the precision of the results presented can be 
inferred from the comparison between L = 95, M = 500, 
and L = 96, M = 300. Clearly, the main features are 
present in both cases but it can be seen that the ones for 
L = 95, M ~ 500 are more consistent as n is reduced. 
For the impurity U' = 8, t' — 0.5, the Kondo peak at 
half-filling remains at the Fermi level upon reducing the 
density down to n = 0.5. This is consistent with the 
relative large conductance of this case. For the impurity 
U' = 8, t' = 1, the peak shifts slightly away from res- 
onance as the filling is reduced, going into the electron 
part of the spectrum (w — /i < 0). Also in this case, 
the behavior of p{0) is correlated with the conductance 
weakly varying with the density. For these two impuri- 
ties, the Coulomb or holon peaks fall outside of the fre- 
quency range adopted in the plot. In the case of U' — 4, 
t' = 0.5, the Coulomb peak at w « —U/2 shows up as 
n decreases. The overall behavior of the Kondo peak is 
similar to that of the first impurity and this is related to 
the small Jeff of these cases. On the other hand, for the 
case U' — 2, t' = 0.5, the chemical potential shifts below 
the Coulomb peak at cj « — C//2 (the Coulomb peaks are 
those located at w — /i ^ ±1 at half-filling). This behav- 
ior is clearly correlated with the important suppression of 
the conductance for this impurity observed in Fig. |2Ia). 
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FIG. 3: (Color online) Evolution of the LDOS at the impu- 
rity site as the electron density is reduced from half-filling to 
quarter-filling. Titles on each panel show the corresponding 
impurity parameters. For each panel, the fillings are n = 1.0, 
0.875, 0.75, 0.625 and 0.5 from top to bottom, L = 95, 
Al — 500. Dashed lines in the top left panel correspond to 
L = 96, M = 300 



For the U' ^ 8, t' = 1, and U' = 2, i' = 0.5 impurities, 
the observed splitting of the central peak at n = 1 is 
just a finite size effect; the split between the two central 
peaks, which correspond to electron creation {uj — p, > 0) 
and annihilation {lu — p < 0), close when L — > cxd. Notice 
also the different amplitudes of the spectral weight at the 
Fermi level, roughly proportional to Jef / . 



IV. INTERACTING LEADS, 
QUARTER-FILLING 

After examining the effect of varying the filling in a sys- 
tem with an Anderson impurity in a tight-binding chain, 
we now turn to the central issue of this paper, that is 
the effect of the presence of a Hubbard on-site potential 
on the chain. Then, we consider the full Hamiltonian 
defined in (Q, at quarter-filling, and taking C/ = 1, 2, 4 
and 8. In the pure Hubbard chain, the system behaves 
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FIG. 4; (Color online) (a),(b) Electron occupancy and (c),(d) 
(Sf ) at the impurity site as a function of the inverse of the 
chain length L for U = Q, 1, 2, 4, and 8, from bottom to top. 
U' = 8, t' = 0.5 (left panel), U' = 2, t' = 0.5 (right panel), 
n = 0.5. 
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FIG. 5: (Color online) (a) Electron occupancy and (b) (Sf ) at 
the impurity site as a function of the interaction on the leads 
U for U' = 8, t' = 0.5 (circles), U' = 8, t' = 1.0 (squares), 
U' = 4, t' = 0.5 (diamonds), and U' = 2, t' = 0.5 (triangles), 
i = 96, n = 0.5. 



as a Luttinger liquid, and hence, from the point of view 
of transport properties is a metal. 

We start by looking at finite size behavior of the elec- 
tron occupancy and (5^) at the impurity. Results for the 
impurities U' = 8, i' = 0.5 and U' = 2, t' ^ 0.5 as a 
function of 1/L arc shown in Fig.^ For each value of U 
a weak variation can be observed in (no). For the case of 
U' = 8, t' = 0.5 there is a somewhat larger variation of 
i^z o) which becomes weaker as U is increased. 

The fact that as U increases, both (no) and {S^^q) in- 
crease, is a general behavior observed for all the values of 
U' and t' that we examined, as can be seen in Fig.[5| Note 
that for U' = 8,t' = 1 in Fig.|S{a) the system reaches the 
limit of the pure Hubbard model at U = 8, with (no) = 1, 
while the other cases represent different kind of impuri- 
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FIG. 6: (Color online) (a) Integral of the spin-spin corre- 
lations as a function of the distance from the impurity site 
(Eq.0), U' = 8, t' = 0.5, and (7 = 0, 1, 2, 4, and 8, from 
top to bottom; L = 96, n = 0.5. (b) Length of the compen- 
sation cloud as a function of U for the four set of interactions 
at the impurity, [/' = 8, = 0.5, (circles), U' = 8, t' ^ 1.0 
(squares), U' = 4, t' = 0.5 (diamonds), and U' = 2, t' = 0.5 
(triangles). L = 64 (open symbols), L = 96 (filled symbols), 
n = 0.5. 



ties, with densities not necessarily equal to 1. This forces 
the crossing of the two U' = 8 curves found near U = 5. 
For the cases U' = 2 and U' = 4, for U > U' we can 
observe a change in the growing behavior of both (no) 
and (5^0) curves. In other words, a Hubbard repulsion 
on the chains favors a more defined magnetic character of 
the impurity. This is particularly apparent for the impu- 
rity with parameters U' — 2, t' = 0.5, which most likely 
is in a mixed valence state at C/ = 0, as discussed in the 
previous Section. Again these results are consistent with 
those obtained with SSE within error bars. One may 
conjecture that as U increases, the Akp electron corre- 
lations are enhanced forcing an increased occupation at 
the impurity. The enhancement of (uq) and {S'^ q) with U 
occurs also at density n = 0.75 at least for this impurity 
parameters. 

Let us examine now the behavior of the spin compensa- 
tion cloud when U is increased. In Fig.EJa), the integral 
of the spin-spin correlations is plotted as a function of the 
distance from the impurity. The typical 2kp oscillations 
which at n = 0.5 has a period equal to 4, are clearly seen. 
S{r) increases (in absolute value) as U increases, imply- 
ing that the impurity spin becomes screened at shorter 
distances. This is precisely the behavior of ^, as it can 
be seen in Fig. EJb) for the four types of impurities con- 
sidered. For U' = 8, t' = 0.5 and U' = 4, t' = 0.5 
the reduction of ^ with U is clear in spite of the fact 
that the dependence of ^ with the system size suggests 
that the bulk limit has not been reached. On the other 
hand, for U' = 8, i' = 1.0 and for U' = 2, t' = 0.5, re- 
sults for L = 64 and 96 are very similar suggesting that 
the compensation clouds observed are those of the bulk 
limit. This behavior of the spin compensation cloud is 
correlated with the increase of (^^ 0) shown in Fig.[S{b). 
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FIG. 7: Evolution of the LDOS at the impurity site as U in 
the leads is increased. Titles on each panel show the corre- 
sponding impurity parameters. For each panel, U — 0.0, 1.0, 
2.0, 4.0 and 8.0 from top to bottom. L = 95, n = 0.5. 



The relative values of ^ for the various impurities fol- 
lows in first approximation the values of the Kondo cou- 
plings, as was discussed with respect to Fig. EJb). The 
overall reduction of ^ with U is consistent with previous 
calculations for the two-impurity Kondo coupling^"^ and 
expected on general grounds due to the increase of the 
Kondo temperature with correlations.^^ 

The LDOS at the impurity presents also interesting 
features as U is increased. Let us first consider the im- 
purity U' ^ 8, t' = 0.5 (top left panel Fig. [7|). In this 
case the peak retains its resonant character as U increases 
(only for f7 = 8 it appears at a small negative frequency) , 
and its amplitude is slightly reduced. For the second type 
of impurity, U' = 8, <' = 1.0 (top right panel Fig. [TI, 
there is a pronounced transfer of spectral weight to the 
occupied part of the spectrum, and also p{uj — = 0) is 
much reduced. For the impurity U' = A, t' — 0.5 (bot- 
tom left panel), the remains of the Kondo peak stays at 
the chemical potential, which moves towards positive w. 
The same effect of U is observed for U' = 2, t' = 0.5 
(bottom right panel Fig. [7||. In this case, the peak lo- 




FIG. 8: (Color online) J{t)/Av for different values of U, (a) 
U' = 8, t' = 1.0, (b) U' = 4,t' = 0.5, L = 96, n = 0.5, AV^ = 
0.01, At — 0.1. The horizontal lines show the time interval 
over which the average of the current is taken to compute 
the conductance, (c) Relative conductance (see text) as a 
function of U for different values of the interactions {U',t'). 
The values of Kp at n = 0.5, extracted from Ref.3^ are shown 
with stars. 



cated at w — /Lt > for J7 = 0, crosses the Fermi level at 
[/ ~ 4, and finally appears at a; — < for J7 = 8, as 
the spectral weight is transfered from the unoccupied to 
the occupied part of the spectrum. Taking into account 
the strong enhancement of (5^ q) for these impurity pa- 
rameters shown in Fig|SJb), one may conjecture that the 
states close to the Fermi level that are shifted by in- 
creasing U correspond to magnetic excitations becoming 
increasingly occupied. 

Our most relevant results regarding the transport 
properties in a Hubbard chain in the presence of an An- 
derson impurity are shown in Fig. |S1 In Fig. |Sfa,b), we 
show the average current J(t) on the links connected to 
the impurity divided by Ay on a L = 96 chain (results 
for L — 95 are virtually identical) for two sets of im- 
purity interactions. As discussed in Section^ J{t) has 
an oscillatory behavior which is typical of these time- 
evolving systemsi^S^ The period of these oscillations for 
the cases shown in Fig. |S1 is approximately 800At. In 
this figure, we only show J{t) up to i = 250At because 
after that the time-evolution becomes unstable. In some 
cases, (Fig.lHJa)), these wiggles make the conductance to 
exceed the value 2e^//i indicating that the time-evolution 
has already become unstable. 



8 



-0.02 



-0.04 



-0.06 



-0.08 



-0.1, 



"T" 



"T" 



"T" 



(a) 



o-oT=0.1 
H-H T=0.075 

T=0.05 
^!"iT=0.025 
< T=0.0125 
• T=0 




40 50 60 70 
r 



_l_ 



_l_ 



80 90 100 




0.06 



FIG. 9: (Color online) (a) Evolution of the spin compensation 
cloud as a function of distance from the impurity at various 
temperatures obtained with SSE, U' — 2, t' — 0.5, L = 96. 
(b) Length of the spin compensation cloud as a function of 
temperature U' = 2, t' = 0.5, L = 96, n = 0.5, for U = 0, 
and 1. The T = values were obtained with DMRG. Error 
bars are shown for T = 0.0125. Temperature in units of t. 



From J{t)/AV, we have computed the conductance ac- 
cording to Eq. (0). Typical time intervals or "plateaus" 
over which the time average was taken are shown in 
Fig. |H^a,b). The values of the conductance G — G{U), 
for a given impurity, relative to the value for noninter- 
acting leads. Go = G{U = 0), are shown in Fig. IHc). 
Although there is some arbitrariness in defining these 
plateaus, a qualitatively clear trend emerges. For the im- 
purities ([/' = 8, t' = 1.0) and {U' ^ A, t' ^ 0.5) G/Go 
decreases as U increases. In these cases, the suppression 
is consistent with the reduction of spectral weight at the 
Fermi level in the impurity LDOS, shown in Fig. [T] In 
these cases, where a Kondo resonance subsists for finite 
U , and particularly for [V — 8, t' — 1.0), the suppres- 
sion of the conductance is consistent with the predicted 
relationi^ G ~ Kp, as shown in Fig. |HIc). On the other 
hand, for the case of the impurity parameters U' = 2, 
t' = 0.5, G/Gq follows a more complex behavior. From 
the non-interacting limit to approximately U = 1, G/Gq 
increases and then as U increases further, it starts to 
decrease. This non-monotonous behavior again follows 
the one of the spectral weight at the Fermi level in the 
impurity LDOS (Fig. [71 bottom right panel). 

Finally, we report a finite temperature study of the 
spin compensation cloud. We chose the impurity defined 
by the couplings U' = 2, t' — 0.5. In this case, since Jeff 
is large (as long as the relation Jeff = At''^/U' is appli- 
cable to these parameters), we expect a relatively small 
length ^ of the spin compensation cloud. In addition, the 
expression for the Kondo temperature for non-interacting 
leads {U = 0), Tk = t'^U'/{2t) exp {--kW / {8t'^)), al- 
though not much precise for these impurity parameters, 
gives Tk ~ 0.02 which is reasonably accessible by SSE 
simulations. In Fig. IHIa) the spin compensation cloud 
S{r) is shown at several temperatures for the case of non- 
interacting leads, [/ = 0. The temperature is in units of t, 



ks = I- It is remarkable the convergence of finite temper- 
ature results obtained with SSE to the zero temperature 
ones obtained with DMRG. Notice that both techniques 
are in principle exact but SSE results correspond to PBC 
systems while DMRG ones to OBC chains. 

The length ^ is shown in Fig. EJb) . Errors were com- 
puted from the deviation of results obtained in four in- 
dependent runs. For T > T*, T* « 0.375, ^ is equal to 
half the chain length, implying that the spin compensa- 
tion cloud exceeds the limits of the chain. For T < T*, 
^ starts to be smaller than half the chain length indicat- 
ing that now the cloud is fitting in the chain. That is to 
say, for T < T* the system manages to compensate the 
spin of the impurity but this is not possible for T > T*. 
If the Kondo temperature is thought as the crossover 
temperature around which the impurity spin becomes 
compensated then T* would be a measure of the Kondo 
temperature in finite systems, that is, T* = T*{L). As 
mentioned in Section m there is a certain degree of ar- 
bitrariness in the estimation of ^ because of the cutoff 
X but what is relevant is to determine how T* behaves 
when a parameter of the system is varied. For the case 
of [/ = 1, T* reduces its value to w 0.25. This seems a 
trend as the correlations on the leads are increased but 
results for larger values of U have larger error bars and 
so they are not precise enough to confirm this trend. 



V. CONCLUSIONS 

In summary, we have studied the single-impurity 
Anderson-Hubbard model on finite chains with numerical 
techniques. In the first place, we analyzed how the reduc- 
tion of the electron density from half-filling to quarter- 
filling affects the Kondo resonance for non-interacting 
leads. Physical quantities defining the magnetic char- 
acter of the impurity show in general a weak behavior 
with the electron density except when the impurity is, at 
half-filling, close to a mixed valence regime. In this case 
there is a steep lost of magnetic character by reducing 
the density. This behavior is refiected also in the im- 
purity LDOS and in the conductance in linear response. 
Next, we turn to the main topic of this work which is the 
study of the effects of a Hubbard interaction on the leads 
at quarter-filling. A very clear and interesting result is 
that U drives the impurity into a more defined Kondo 
regime (Fig.O although accompanied in most cases by a 
reduction of the spectral weight of the impurity LDOS. 
This reduction is however not general since for an im- 
purity close to the mixed valence regime we observed a 
nonmonotonic behavior. As we mentioned in the Intro- 
duction, there is a number of strongly correlated materi- 
als which can be or have already been implemented- in 
nanoscopic devices and where some of the present results 
for the impurity LDOS could be contrasted. The con- 
ductance, computed for a small finite bias applied to the 
leads, follows the behavior of the impurity LDOS. For the 
impurities in a well-defined Kondo regime, the conduc- 
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tance is suppressed by U in finite chains, in agreement 
with what is expected for an impurity in a Luttinger hq- 
uid. Again, for the impurity close to a mixed valence 
regime, G increases for small U but eventually starts to 
decrease for large values of U. This nonmonotonic behav- 
ior indicates a complex interplay between the interactions 
on the impurity and the ones on the leads (in addition, 
of course, to the impurity-leads interactions). 

Finally, we exploit one of the important advantages of 
the real space methods we used, that is, the possibility of 
getting a measure of the spin screening or compensation 
of the impurity. The "spin compensation cloud" gives an 
insight on the internal working of the Kondo effect and 
it does not depend just on properties at the impurity site 
but also on U on the leads. As shown in Fig. finite- 
size effects are more important for small values if U. The 
proposed determination of the Kondo temperature in fi- 
nite systems through the spin compensation cloud, has 



to be thoroughly explored as a function of density and 
coupling parameters. We emphasize that this numeri- 
cal estimation of T* is completely general in the sense 
that it is not limited to ID Anderson-Hubbard model 
but could be used for any Hamiltonian for the leads, for 
example the Kondo lattice model, and also not limited to 
ID leads, although beyond ID the "minus sign problem" 
makes difficult the application of QMC methods. More- 
over, we believe that the dependence of T* with the finite 
size of the system could be experimentally measured and 
compared with theoretical results. 
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